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In a more general electroweak theory, there could be Higgs particles that are odd under 
CP. When such particles decay via vector bosons to two fermion-antifermion pairs, the 
momenta of those will be correlated in a way which is determined by the CP of the particle. 
Similarly, in the Bjorken process correlations among momenta of the initial electron and 
final-state fermions are sensitive to the CP quantum number. Monte Carlo data on the 
expected efficiency demonstrate that it should be possible to verify the scalar character of 
an intermediate-mass Standard Model Higgs boson after three years of data taking at a 
future linear collider. This is most likely not possible at LEP2. Signals of possible presence 
of CP violation in the Higgs sector are briefly discussed. 
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1 Introduction 



The only known CP violation is in the .fT-meson system, where so far five CP violating 
observables have been studied. These are the three "classical" ones pQ 
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and two more recently observed ones [2HH] 
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One hopes to also study CP violation in the B meson system, at the planned SLAC B 
factory and dedicated future detectors at other facilities like HERA and the LHC The 
above observables (which are very small) are often expressed in terms of e and e', e.g., 

?7 + _ ~ e + e 
r] 00 ~ e - 2e' 
5i ~ 2 Re e (3) 

with |e| = (2.27 ± 0.02) ■ 10" 3 and |e'| < e. 

It has been known since 1967 j3] that CP violation is a necessary prerequisite for 
creating a matter-antimatter asymmetry, as the Universe appears to have. But while 
the CKM matrix of the standard model naturally accomodates CP violation through the 
complex phase factors in charged-current couplings of the quarks 0, it appears that CP 
violation through the CKM matrix may not be sufficient to explain the baryon asymmetry 
of the Universe j^j, so it is important to explore alternative sources. Many years ago 
Weinberg suggested that CP might be violated in the Higgs sector, and presented an 
explicit model where this is the case [Zj. 

We shall here explore the possibility of CP violation in the Higgs sector, following a 
purely phenomenological approach. Of course, we assume a non-minimal Higgs sector. 
Assuming one finds some Higgs-particle candidate, there are two related questions of fun- 
damental interest: (z) Is the Higgs candidate even or odd under CP1 and (ii) Is CP 
violated in the Higgs sector? 

Two physical situations allow us to address these questions in a phenomenological way: 

H^V 1 V 2 ^{f 1 f 2 ){hU) (LHC) (4) 

e+e- - ZH - {fJi)H (NLC) (5) 
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In several ways, they are quite similar. 

We adopt the following effective Lagrangian describing the HVV-vertex: 

C HVV = 2 ■ 2 1 ' A ^[m 2 v V% H + \ V V,„V pa H] (6) 

where = d^Vv — dyV^. The first term describes the standard-model HVV coupling 
(even in CP), whereas the second term describes the coupling of a CP-odd scalar to two 
vector bosons. (In the MSSM jH], this latter coupling is absent at tree level.) Armed with 
this Lagrangian, we address the above two questions in the following two sections. 

2 CP Eigenstates 

Let us first consider Higgs decay, as it could presumably be studied at the LHC: 

H - v x V 2 - [/i(gi)/ 2 (?2)][/ 3 (?3)/ 4 (g4)] (7) 

where ViV 2 denote W + W~ or ZZ . The differential decay rate can be written as |§HTT] 



d 8 r, = q 



Xi + sin(2xi) sin(2x2)^ dLips(m 2 ; q u q 2 , q 3 , qi), i = H,A (8) 



where H denotes a CP-even Higgs particle and A a CP-odd one. Furthermore, 

Xu = (qi ■ q s )(q2 ■ Qi) + (qi ■ qd(l2 ■ qs) 
Yr = (qi ■ q 3 )(q 2 • q*) - (qi ■ g 4 )(<?2 • qa) 
X A = -2[( qi -q 2 )(q 3 -q 4 )] 2 

-2[(ft • gs)(?2 • qi) - (qi ■ qi)(q2 ■ gs)] 2 

+ (qi ■ q2)(q3 ■ q*){[(qi ■ ?a) + (?2 ■ qi)} 2 + [(qi • qi) + (?2 • gs)] 2 } 

Y A = (qi ■ <?2)(<?3 ■ qdiiqi - <h) ■ + qi)}[(qs - qd ■ (qi + 92)] (9) 

These decay rates depend on the relative azimuthal orientation of the planes formed 
by the decay products, in a way which reveals the CP quantum number of the decaying 
particle. It is natural to define an azimuthal angle by 

, (qi X q 2 ) ■ (q 3 X q 4 ) 
COS0=— r: j— (10) 

|qi x q 2 ||q3 x q 4 | 

and consider the corresponding decay distributions. (This is completely analogous to the 
classic determination of the tt° parity [EJ-) Integrating over the way the energy is shared 
within each pair, over the invariant masses of the two pairs and the polar angle distribution, 
we are left with a dependence on the azimuthal angle <fi only. The term Ya does not 
contribute to the resulting decay distributions, which can be expressed in the following 
compact way jOHHj 

— = 1 — a{m) sin(2%i) sin(2x 2 ) cos0 + (3{m) cos20 

r H d0 

^ = l-icos20 (11) 

r A d(p 4 
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where the sin(2x) factors are given by ratios between the axial and vector couplings [TU] . 
The coefficients a(m) and (3{m) depend on whether the decay takes place via W + W~ 
or via ZZ, and specifically on the ratio of the vector boson mass to the Higgs mass m. 
Representative distributions ^0] are shown in fig. 1 (left part). Of course, jet identification 
is needed for this kind of analysis. 
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Figure 1. Left part: Characteristic angular distributions of the planes defined by two Dalitz pairs 
for CP-even Higgs particles decaying via two W's or two Z's, compared with the corresponding 
distribution for a CP-odd Higgs particle (denoted A). Different Higgs masses are considered in 
the CP-even case. Right part: Characteristic energy-weighted distributions for CP-even Higgs 
particles decaying via two W's or two Z's, compared with the corresponding distribution for a 
CP-odd Higgs particle (A). (Ref. \10)j had a sign error in these distributions.) 

In order to project out the YA-term of eq. 0, we multiply (JBJ by the energy differences 
(oji — oj2){u)3 — before integrating over energies. The distribution then takes the form 

DDI 



2tt dr 



1 - 



k (m) 



COS( 



'121 



sin(2xi) sin(2x 2 ) 

in the CP-even case, whereas it is constant, nim) = 0, in the CP-odd case. Illustrative 
distributions are given in fig. 1 (right part). This energy- weighted angular distribution 
provides an independent determination of the CP of the Higgs candidate, and is in fact 
more sensitive than the unweighted angular distributions, especially for heavy Higgs bosons. 
Consider next the Bjorken process [13] 

e (Pi) e + (pa) -+Z(Q)h (q 3 ) - / ( 9l ) / (q 2 )h (q 3 ) (13) 
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Let the momenta of the two final-state fermions and the initial electron (in the overall cm. 
frame) define two planes, and denote by the angle between those two planes; i.e. 



cos « 



vPi x Q) • (qi x Q2) 



|Pi x Qllqi x q 2 | 
We shall discuss the angular distribution of the cross section a, 

2tt da 
a d(f> 



(14) 



(15) 



both in the case of CP-even and CP-odd Higgs bosons. These distributions take the form 

HE EU El EE] 



2n dan 
a H d(f) 
2n ddA 
a A d(p 



1 + a(s, m) cos <f> + /3(s, m) cos 20, 



1 COS : 
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(16) 
(17) 



with typical cases shown in fig. El There is a clear difference between CP = 1 and CP = —1 
for a wide range of energies and Higgs masses. 
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Figure 2. Angular distributions of the planes defined by incoming e~ and final-state fermions 
for a CP-even Higgs particle (solid) compared with the corresponding distribution for a CP-odd 
one (dashed). Different energies and masses are considered in the CP-even case. We assume 
i/i = 200 and 500 GeV at LEP2 and NLC, respectively. 



Experimentally, however, one faces the challenge of contrasting two angular distribu- 
tions with a restricted number of events and allowing also for background. We shall here 
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focus on the intermediate Higgs mass range; more specifically, we consider m < 140 GeV 
where the Higgs decays dominantly to bb. The main background will then stem from 
e + e~ — > ZZ and also e + e~ — > Zj, 77. The cleanest channel for isolating the Higgs signal 
from the background is provided by the and e + e~ decay modes of the Z boson. 

Let us next limit consideration to the energy range \fs = 300 — 500 GeV, as appro- 
priate for a linear collider ^3], henceforth denoted NLC. We impose reasonable cuts and 
constraints as described in [TB]; e.g. \m IM + IM - —m z \ < 6 GeV and | cos 0^1 < 0.6. The signal 
for e + e~ — > ZH — > fi + [i~bb will then be larger than the background e + e~ — > ZZ — > fi + {i~bb 
by an order of magnitude. In the following we shall thus neglect the background in the 
discussion of fTB]) versus (fTTj) . With cx(e + e" — > ZH) ~ 200 fb and an integrated lu- 
minosity of 20 fb" 1 a year [IB] , about 4000 Higgs particles will be produced per year, 
in this intermediate mass range. However, following |TB| we have only ~ 30 signal events 
e + e~ — > ZP — > ji + ji~bb left per year for e.g. a NLC operating at y/s = 300 GeV and a Higgs 
particle of mass m = 125 GeV. In the case e + e~ — > ZH —>■ e + e~bb we also have a t-channel 
background contribution from the ZZ fusion process e + e~ — > e + e~(ZZ) — > e + e~H. This 
contribution may be neglected at LEP energies, but it is comparable to the s-channel con- 
tribution at higher energies. However, this contribution can be suppressed by imposing a 
cut on the invariant mass of the final-state electrons, e.g. \m e + e - — mz\ < 6 GeV. Hence, 
we can effectively treat the electrons on the same footing as the muons, thereby obtaining 
a doubling of the event rate. 

In order to demonstrate the potential of the NLC for determining the CP of the Higgs 
particle, we show in fig. 3 the result of a Monte Carlo simulation. For this purpose we 
have used PYTHIA [17] . suitably modified to allow for the CP-odd case. The statistics 
correspond to 3 years of runnings using both the and e + e~ decay modes of the Z 

boson. This yields about 200 events in these channels, assuming an integrated luminosity 
of 20 fb -1 per year. 

For -y/s = 300 GeV and m# = 125 GeV, the 'bare' prediction for f3 is 0.12. The cut 
on polar angle, | cos 6*| < 0.6, increases it slightly to 0.14. Similarly, the ' — 1/4' changes 
to —0.39. Thus, the cut makes it easier to discriminate between the CP-even distribution 
and the CP-odd one [H]. From fig. 3 we see that the individual angular Monte Carlo 
distributions are consistent with the predictions, showing that a three-year data sample is 
large enough to reproduce the azimuthal distributions. In the Standard-Model case the fit 
gives 0.92±0.07 and 0.2±0.1 for the predictions 1.00 and 0.14, respectively, with \ 2 — 1-0. 
In the CP-odd case the fit gives 0.94 ± 0.07 and —0.4 ±0.1 for the predictions 1.00 and 
—0.39, respectively, with x 2 — 0.7. More importantly, since the cos 20 terms are more 
than 4 standard deviations away from each other, a data sample of this size is sufficient to 
verify the scalar nature of the Standard- Model Higgs. Using likelihood ratios, as described 
in [THj, for choosing between the two hypotheses of CP even and CP odd, we find that 
less than 3 years of running suffices. 

An alternative test has recently been suggested by Arens et. al. ^2] in the context of 
Higgs decaying via vector bosons to four fermions, where one studies the energy spectrum 

■i-The event rate is based on the Standard Model, and could be different for a non-standard Higgs sector. 
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Figure 3. Monte Carlo data displaying the angular distribution of events e + e~ — > ZH — > l + l~bb, 
I = fj,, e for a Standard-Model Higgs versus a CP-odd one. We have taken = 300 GeV, 
m = 125 GeV, and an angular cut \ cos6\ < 0.6. 



of one of the final-state fermions. Applying this idea to the Bjorken process one would 
study the energy distribution of an outgoing fermion, e.g. u\~ or e~. Introducing the scaled 
lepton energy, x = AEi- / y/s, I = /i, e, we shall consider the energy distribution of the cross 
section with respect to this final-fermion energy, 

i £ (18) 
a ax 

both in the case of CP-even and CP-odd Higgs bosons. We are using the narrow-width 
approximation and the range of x is given by x_ < x < x + , with sx± = s + m| — m 2 ± \f\. 
Here the distributions are given as second-degree polynominals in x, and, as shown in [Ti] . 
the coefficients have a non-trivial dependence on the cm. energy and the Higgs mass, also 
for the CP-odd case. A representative set of energy distributions is given in fig. 4 for the 
case e + e~ —* ^ nrh for both LEP2 and NLC energies. Again, there is a clear difference 
between the CP-even and the CP-odd cases. 

A Monte-Carlo simulation for the energy distribution eq. ()18|) confirms this statement 
|14j . An analysis of likelihood ratios demonstrates that less than 3 years of running is 
sufficient if we require a discrimination by four standard deviations, but more events seem 
to be required than in the case of angular distributions. 



3 CP Violation 

If we allow for both the Standard- Model and the CP-odd term in the Higgs- vector coupling 
(jnj), then there will be CP violation in Higgs decay of the kind (@J) as well as in the Bjorken 
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Figure 4. Characteristic distributions for the scaled energy, eq. ifT%|) . of the I , I = fi,e in the 
Bjorken process e + e~ — > l + l~h. Different energies and masses are considered. 



process [UJ HH [TBI EDI 121] • The decay rate can be written as 



d 8 r = Vl — N^D 
m 



X + \ Y +(^T^ Z 



7717 



in 



v 



dLips(m 2 ;g 1; 42,43,44) (19) 



where the Y term describes CP violation. 

Integrating over the way in which energy is shared in each fermion-antifermion pair, as 
well as over polar angles, we can write the distribution in a compact form as [21] 



d 3 r 



v " N 1 N 2 ^ y /X(m^s 1 ,s 2 )D(s 1 ,s 2 ) 



d(/>dsids 2 72(4tt) 6 

A (m 2 , s u s 2 ) + 4sis 2 (1 + 2p 2 ) + 2siS 2 p 2 cos 2(0 - 5) 



x 



3^ 



- sin 2%i sin 2%2 ( — ) y/s\s 2 (m - s 1 - s 2 ) p cos(0 - 5) 



with a modulation function p 



p = ^1 + 7] 2 X (m 2 , si, s 2 ) / (4m£) 
and a "tilt" given by the amount of CP violation as 

77(si, s 2 )^JX(m 2 , si, s 2 ) 



5 = arctan ■ 



2m 2 , 



tt/2 < (5 < tt/2 



(20) 



(21^ 



(22) 



and X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) the Kallen function. 
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An inclusive distribution is obtained if we integrate also over the invariant masses of 
the two pairs. It has the form 

— — = 1 - aim) p cos(0 - 5) + f3(m) p 2 cos 2 (0 - 6) (23) 

r d</> 

The tilt 5 will depend on the mass of the Higgs particle, and is shown for some represen- 
tative values of rj in fig. 5 (left part). If the Higgs is somewhat heavy, and rj not too small, 
this tilt can be significant. Again, energy- weighted correlations provide an independent 
way of measuring this kind of CP violation [21]. 

The same kind of analysis may be applied to the Bjorken process. There is a relative 
shift in the angular distribution of the two planes, defined by eq. (JHJ), due to CP violation. 
This rotation vanishes at the threshold for producing a real vector boson and, even for a 
fixed value of rj, grows with energy. The distribution takes the compact form 

— ^ = 1 + a(s, m) p cos(0 + 5) + /?(s, m) p 2 cos 2 (0 + S) (24) 
a dcp 

Any CP violation would thus show up as a tilt in the azimuthal distribution, by the amount 
5. " 




Figure 5. Left part: The angle 5 (in degrees) for a Higgs particle of mass m, for rj = 1, 
10 _1 , 10~ 2 , and 10 -3 . Right part: Characteristic angular distributions for different amounts 
of CP violation in the Bjorken process, including the CP-even (rj = 0) and CP-odd (rj S> 1) 
eigenstates. We have used rj = 0.1, 0.5, 5 for ^/s = 500 GeV and m = 200 GeV. 

A representative set of angular distributions is given in fig. 5 (right part) for a broad 
range of rj values. We have considered a Higgs boson of m = 200 GeV accompanied by 
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a /i + yU~-pair in the final state, produced at \fs = 500 GeV. We observe that for i] < 0.1 
and r] > 5, the deviations from the CP-even and CP-odd distributions, respectively, are 
small. Experimentally it will be very difficult to disentangle two distributions which differ 
by such a small phase shift. Thus, observation of a small amount of CP violation would 
require a very large amount of data. 

We note that the special cases rj = and rj ^> 1 correspond to the CP even and CP 
odd eigenstates, respectively. Hence, the distribution (j2jj) should be interpreted as being 
intermediate between those for the two eigenstates. 

4 Summary 

We have discussed ways to analyze the CP properties of Higgs candidates through angular 
decay distributions. In particular, we have addressed the problem of estimating the amount 
of data needed in order to distinguish a scalar Higgs from a pseudoscalar one at a future 
linear collider. This is most likely not possible at LEP2 due to much smaller event rates 
and background which is not easily suppressed. However, we have demonstrated that one 
will be able to establish the scalar nature of the Higgs boson at the Next Linear Collider 
from an analysis of angular or energy correlations. This particular study has been carried 
out for the case \fs = 300 GeV, m = 125 GeV. Similar results are expected in other cases 
as long as the background is small. 

In order to establish or rule out specific models, one will also need to compare different 
branching ratios, in particular to fermionic final states. The methods discussed above 
instead deal with quite general properties of the models. 

We are grateful to the Organizers of the Third Tallinn Symposium on Neutrino Physics, 
in particular Professor I. Ots, for creating a stimulating and pleasant atmosphere during 
the meeting. This research has been supported by the Research Council of Norway. 

References 

[1] CP Violation, edited by L. Wolfenstein, North-Holland, 1989. 

[2] E.J. Ramberg et ai, Phys. Rev. Lett. 70 (1993) 2529; J.N. Matthews et ai, Phys. 
Rev. Lett. 75 (1995) 2803. 

[3] CPLEAR Collaboration, R. Adler et ai, CERN PPE/95-134. 

[4] A. D. Sakharov, JETP Lett. 5 (1967) 24. 

[5] M. Kobayashi and M. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 

[6] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Annu. Rev. Nucl. Part. Sci. 43 (1993) 
27. 



10 



[7] S. Weinberg, Phys. Rev. Lett. 37 (1976) 657. 

[8] P. Fayet, Phys. Lett. 69B (1977) 489; 

S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150; 

H.-P. Nilles, Phys. Rep. C110 (1984) 1; 

H. E. Haber and G. L. Kane, Phys. Rep. C117 (1985) 75. 

[9] C. A. Nelson, Phys. Rev. D37 (1988) 1220; 

M. J. Duncan, G. L. Kane and W. W. Repko, Phys. Rev. Lett. 55 (1985) 579; Nucl. 

Phys. B272 (1986) 517; 

B. A. Kniehl, Phys. Lett. B244 (1990) 537. 

[10] A. Skjold and P. Osland, Phys. Lett B311 (1993) 261. 

[11] K. Hagiwara and M.L. Stong, Z. Phys. C 62 (1994) 99; 

M. Kramer, J. Kiihn, M.L. Stong and P.M. Zerwas, Z. Phys. C 64 (1994) 21; 
M.L. Stong, Invited talk, Ringberg Workshop "Perspectives for Electroweak Interac- 
tions in e + e" Collisions", Munich, Feb. 5-8, 1995, hep-ph/9504345, 

[12] C.N. Yang, Phys. Rev. 77 (1950) 242, 722; 

R.H. Dalitz, Proc. Phys. Soc. (London) A64 (1951) 667; 
N. Kroll and W. Wada, Phys. Rev. 98 (1955) 1355. 

[13] J. D. Bjorken, in Proceedings of the 1976 SLAC Summer Institute on Particle Physics, 
ed. M. Zipf (SLAC Report No. 198, 1976) p. 22; 
B.L. Ioffe and V.A. Khoze, Sov. J. Part. Nucl. 9 (1978) 50; 
D.R.T. Jones and S.T. Petcov, Phys. Lett. 84B (1979) 440; 
J. Finjord, Physica Scripta 21 (1980) 143. 

[14] A. Skjold and P. Osland, Nucl. Phys. B 453 (1995) 3. 

[15] B.H. Wiik, in HEP 93, Proc. Int. Europhysics Conference on High Energy Physics, 
Marseille, 22-28 July, 1993 eds. J. Carr and M. Perrottet (Editions Frontieres, Gif- 
sur-Yvette, 1994) p. 739. 

[16] V. Barger, K. Cheung, A. Djouadi, B.A. Kniehl and P.M. Zerwas, 
Phys. Rev. D49 (1994) 79, 

B.A. Kniehl, Proc. Workshop on Physics and Experiments with Linear e + e~ Colliders, 
Waikoloa, Hawaii, April 26-30, 1993 vol. II, p. 625. 

[17] T. Sjostrand, Computer Physics Commun. 82 (1994) 74. 

[18] B.P. Roe, Probability and Statistics in Experimental Physics, ( Springer- Verlag Inc., 
New York, 1992). 

[19] T. Arens, U. D. J. Gieseler and L. M. Sehgal, Phys. Lett. B339 (1994) 127, 
T. Arens and L. M. Sehgal, Z. Phys. C66 (1995) 89. 



11 



[20] J. R. Dell'Aquila and C. A. Nelson, Phys. Rev. D33 (1986) 101; 
B. Grzadkowski and J.F. Gunion, Phys. Lett. B294 (1992) 361; 
D. Chang, W. -Y. Keung and I. Phillips, Phys. Rev. D48 (1993) 3225. 

[21] A. Skjold and P. Osland, Phys. Lett. B329 (1994) 305. 



12 



